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Abstract 

The thermalization rate of a heavy quark is related to its momentum diffusion coefficient. 
Starting from a Kubo relation and using the framework of the heavy quark effective theory, we 
argue that in the large-mass limit the momentum diffusion coefficient can be defined through 
a certain Euclidean correlation function, involving color-electric fields along a Polyakov loop. 
Furthermore, carrying out a perturbative computation, we demonstrate that the spectral 
function corresponding to this correlator is relatively flat at small frequencies. Therefore, 
unlike in the case of several other transport coefficients, for which the narrowness of the 
transport peak makes analytic continuation from Euclidean lattice data susceptible to severe 
systematic uncertainties, it appears that the determination of the heavy quark thermalization 
rate could be relatively well under control. 
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1. Introduction 



One of the very interesting discoveries of the RHIC program at Brookhaven National Labo- 
ratory has been that heavy quarks (particularly the charm quarks) appear to thermalize just 
about as effectively as the light quarks. This has been inferred from measuring the electron 
PT-spectra produced by the decays of the heavy quarks, showing indications of the same type 
of hydrodynamic flow as experienced by the lighter quarks [Ij. 

On the theoretical side, the historical starting point for a QCD-based understanding of the 
behavior of heavy quarks in a thermal environment was the determination of their energy loss 
rate to the leading order in the QCD coupling constant, ctg [2]. The energy loss is directly 
related to a number of other concepts, such as the diffusion and the thermalization rates of 
the heavy quarks [31 [4]. In particular, assuming that the effective value of as is relatively 
small leads to the thermalization rate T ~ a^T'^/M, where T is the temperature and M is 
the heavy quark mass [1]. The comparable thermalization rate for a light quark or gluon is 
r ~ a^T, or at very high energies T ~ a^T s/tJE [5]. Hence heavy quarks with M ^ T are 
expected to thermalize slowly, particularly at weak coupling. 

As already mentioned, the empirical facts appear however to be in conflict with a slow 
thermalization rate [1]. This has lead to a lot of new theoretical ideas, with the hope of 
bringing the theoretical determination of T beyond the leading order in Og. In particular, 
possibilities for a lattice determination were explored [6] ; computations in a strongly coupled 
theory with some similarities with QCD were carried out [71 [8]; phenomenological model 
treatments of bound states were considered [9]; and the first weak-coupling corrections to 
the leading order result were determined [10]. The studies [7]-[l0] showed that there could 
indeed be substantial corrections (of a positive sign) to the leading order result; at the same 
time, the study [6] showed that a direct lattice determination of the heavy quark related 
observables would be extremely hard, because the physics resides in a "transport peak" of 
a certain spectral function, of width Au; ~ F ~ Og <C T, which regime is practically 

impossible to explore with Euclidean techniques. 

The purpose of this paper is to reconsider the prospects for a lattice determination, making 
use of heavy quark effective theory [llj in order to systematically consider the behavior of 
the heavy quarks in the limit M ^ T. Essentially, this allows us to restrict the attention to 
the "numerator" of the thermalization rate, ~ a 2r2, which remains finite in the heavy quark 
limit. In fact, our main goal will be to derive an observable measurable on the lattice which 
has its structure at "large" frequencies, uj ~ a^^T » F, and can be addressed much more 
easily than F itself. 

We note that, in many respects, our study parallels that in ref. [8]. The main differences 
are that we use the imaginary-time formalism rather than the real-time one, in order to make 
contact with the Euclidean spacetime accessible to lattice techniques; and that we try to keep 
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explicit track of 0{as) quantum corrections and renormalization issues. 

The plan of this paper is the following. In Sec. [2] we derive, by going through several 
intermediate steps, the observable alluded to above, consisting of color-electric fields along 
a Polyakov loop. In Sec. [3] we analyze the corresponding spectral function perturbatively, 
demonstrating a relatively flat behavior at small oj < al^^T. In Sec. H] we suggest a lattice 
discretization for the object derived in Sec. El while Sec. [5] offers some conclusions and an 
outlook. 

2. Reduction of the current-current correlator 

In order to proceed with our derivation, we focus on one of the heavy quarks of physical QCD; 
either the charm or the bottom quark. We assume for the moment the use of dimensional 
regularization in order to regulate the theory (though we do not indicate this explicitly). 
Then there is only one large scale in the system, namely the (renormalized) heavy quark 
mass, M, and the task is to account for its effects analytically in the basic observable to be 
defined presently (Eq. ()2.ip ). In Sec. H] we return to the complications emerging in lattice 
regularization. 

2.1. Definitions 

The diffusive motion of heavy quarks within a thermalized medium can be characterized by 
four different quantities, all of which are related to each other (at least in the weak-coupling 
limit). We start by defining the "diffusion coefficient", D, proceeding then to the "relaxation 
rate" or "drag coefficient", rju, and the "momentum diffusion coefficient", n. The fourth 
quantity, the energy loss dE/dx, is historically the first one addressed within QCD pj; yet 
it is not obvious how it could be related to the others on the non-perturbative level, so we 
omit it from the discussion below. 

Among the three quantities, the one that can most directly be defined within quantum field 
theory is the diffusion coefficient D. Consider the spectral function related to the current- 
current correlator, 

(.;) = /"die-*! d^x [J^(t,x), ^^^(0,0)]^ , (2.1) 

where = ip ip; ip is the heavy quark field operator in the Heisenberg picture; (...) = 
Z~^Tr [{...)e~^^] is the thermal expectation value; and (3 = 1/T is the inverse temperature. 
Diffusive motion leads to a pole in the spectral function at a; = —iDh^, where k is the 
momentum (already set to zero in Eq. (j2.ip l. Solving for the pole position and making use 
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of various symmetries leads to the Kubo relation (see, e.g., chapter 6 of ref. [12 



1 

D = — lim > 

3^00 ^^0^ 



i=l 



(2.2) 



Here 'x^^ corresponds to a "susceptibility" related to the conserved charge Jd^x j'''(t,x). 



(2.3) 



For a dilute system of heavy quarks, defines their "number density"0. Note that the 

conserved vector current does not require renormalization, so that the definitions in 
Eqs. (j2.2p and (j2.3p are guaranteed to be ultraviolet finite at any order. 

To define the other quantities, we need to assume that the spectral function around zero 
frequency possesses a narrow transport peak. Due to a heavy quark's large inertia, this is 
certainly true for M sufficiently large, which we assume to be the case from now on. In this 
limit, the spectral function will on general grounds take the form of a Lorentziari, 
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(2.4) 



where Wuv is a frequency scale at which the Lorentzian is overtaken by other types of physical 
processes. 

The other two transport coefficients are then defined from the properties of the transport 
peak. We define the "drag coefficient" r]D to be the width of the Lorentzian, and the (a priori 
mass-dependent) "momentum diffusion coefficient" k^^^ to be -M^j^^ times the coefficient of 
the power-law falloff of its tails. 
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Here Mkin refers to the heavy quark's kinetic mass, to be defined presently (cf. Eq. (j2.7p ). 
Later on we will define a transport coefficient k from the M — > oo limit of k^^^ . 

The physical motivation for the definition in Eq. (j2.5p is as follows. In the dilute limit the 
current couples individually to the heavy quarks; the spectral function [yS) is thus a 
product of their number density times a contribution from one heavy quark. For a single 
quark, J d'^x^j* = represents a non-perturbative measurement of its velocity. Recalling 
Newton's law, AfkindiTVdi is the force acting on the heavy quark; thus k^^^^ is a correlator 



In the non-relativistic limit and at zero chemical potential, Tx"" = ^N^iMT 12-Kfl'^ exp(-/3M); however, 
our basic arguments hold also at a non-zero chemical potential for the heavy quarks, whereby the exponential 
suppression could be removed from T^'"'- 

^Two concrete examples for how such a dependence on uj can arise are reviewed in appendix A. 
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of that force with itself at different times, transformed into frequency space. The factor 
2T/uj relates the spectral function to a time-symmetric correlator, for which this classically 
motivated argumentation applies. Thus Eq. (j2.5p generalizes, in a non-perturbative way, the 
force- force correlator definition of k given in ref. [8]. The condition on w instructs us to 
integrate this force over a time scale, long compared with the medium's correlation time (set 
by t ~ tL>uv)) but short compared with the dynamics of the heavy quark (set by t ~ t]'^). 

The coefficients D, rjD and k^^^ thus defined are related by fiuctuation-dissipation relations, 
which follow from the fact that the area under the transport peak defines the (coarse-grained) 
equal-time mean-squared velocity of a heavy quark, 

(V ) = ^ y_ ^ L ^^G„v(^) . (2.6) 

Here we have introduced a cutoff function GuYiio) designed to isolate the transport peak from 
other types of physics, for instance Guv('^) = ^(^uv — |w|). In the time domain we are thus 
averaging over a time scale tuv ^ "^uv- Such a time averaging is mandatory to make (v^) 
finite and well-defined, since an instantaneous measurement of the heavy quark's velocity 
would induce it to radiate (or absorb) large amounts of energy, thereby changing its state. 
We emphasize that, were there no sharp zero-frequency peak in the spectral function, there 
would be no unambiguous notion of a heavy quark's (coarse-grained) mean squared velocity. 

Motivated by the standard non-relativistic thermodynamic result, we can now define a 
kinetic mass via 

(v') ^ 3-^ . (2.7) 
Eqs. (|2.4p - (j2.6p then yield the fluctuation-dissipation, or Einstein, relations: 

2r2 k(a^) , , 

both of which involve 0{rj£)/uj^y) relative uncertainties. Note that r]£) ~ 1/Mkin in the large 
mass limit, assuming that k^^^) contains no (power-like) dependence on M^in, justifying the 
narrow peak assumption. 

Thermodynamic considerations relate the kinetic mass defined in Eq. (|2.7p to the standard 
notion. Namely, thanks to the slow dynamics of a heavy quark, one can (approximately) 
define a free energy F{v) as a function of its velocity (time-averaged over a period ~ tuv); 
expanding it as -F(v) = Mrest + -^kinV^/2 + ©(v^) at small v and taking a thermodynamic 
average should reproduce Eq. (j2.7p . 

The only approximations we have made so far concern the assumption of a narrow transport 
peak. Parametrically, in weakly coupled QCD [3], pY{io)/{x^^uj) has a peak value D ~ l/g'^T, 
where = 47ras; a width rjD ~ g'^T'^/M; and a perturbative ultraviolet contribution which 
will start to depart from the Lorentzian tail at the scale ~ gT (see Sec. [3]). Thus 
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errors are of order rj£)/gT ~ g'^T/M. In strongly coupled multicolor {Nc oo) Af = 4 Super- 
Yang- Mills theory [7], with a 't Hooft coupling A = g'^Nc, the width of the transport peak is 
Vd ^ ^/XT^/Ad, and the continuum takes over at Wuv ~ T; thus ambiguities are suppressed 
by y/XT/M. 

Expecting the force- force correlator k*^^^) to actually be mass-independent at large Mkin, 
as will be verified a posteriori, we are finally led to take the Mkin oo limit of Eq. (|2.5p . 
inside of which it is essential to retain lo small but non-zero: 



— lim uj"^ 
1=1 



kin 



lim 



00 



dte-yd3x^-{j^(t,x),J^(0,0)} 



(2.9) 



The factor 2T/uj has been accounted for by replacing the spectral function by a time- 
symmetric correlator. Eq. ()2.9p will be the starting point for the further steps to be taken. 



2.2. Heavy quark limit 

Starting from the definition in Eq. (|2.9|) . our next goal is to carry out the limit M — > cxd. As 
a first step we note that, making use of time translational invariance and carrying out partial 
integrations, the definition in Eq. (|2.9|) can be rephrased as 



K = —y lim 

3 ^c^^o 
1=1 



lim — l^'" 

A/^oo X 



00 



dt 



dt' 



(2.10) 



In order to evaluate the time derivatives here, let us rewrite the QCD Lagrangian, £qcd = 
^5(17^1)^4 — M)7p + flight, after a Foldy-Wouthuysen transformation [13]: expanding in 1/M 
and dropping total derivatives, this yields 



/:qcd = e^(iDo-M+'-^^^^^p^y + c^^[iD,+ 



M 



C2 + eg q- • gB 
2M 



'^{0^a-gBcj)-^^a-gBi 



(2.11) 



where Di = di — igAi, gBi = ^€ijk[Dj,Dk], gEi = i[Do,Di], and 6,(j) are two-component 
spinors. The mass M is the pole masqj. 



M = m{fl)\ 1 + 



(47r)2 



In 



^ _ + - 
m?{^) 3 



+ 0{g^ 



(2.12) 



where m(/i) is the MS mass. In regularization schemes respecting Lorentz invariance, the 
coefficient C2 must equal unity (because the combination needed for solving for the pole mass 



^ This is true in schemes producing no additive mass renormalization, such as dimensional regularization. 
There is no multiplicative renormalization to M in Eq. (|2.1H) either, because M could be shifted to zero by 



the field redefinitions 



and would then remain zero quantum mechanically. 
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is ~ Pq — — M^), and we assume this to be the case in the following. The matching 
coefficients cb,ce equal unity at leading order but have quantum corrections; these are not 
needed in the present study. Note that the linearly appearing "rest mass" M is normally 
shifted away (or rather replaced with 0"^); however, we prefer to keep it explicit for the 
moment, because the shifts needed are non-trivial at a non-zero temperature, where the 
Euclidean time extent is finite. 

Setting C2 = 1 in Eq. ()2.1ip . we can read off the conserved Noether current and the 
Hamiltonian in the heavy quark- mass limit: 



J' 
J' 

H 



2M I 
d^x 



(2.13) 
(2.14) 

(2.15) 



Here we treat the fermionic fields as operators but the gauge fields as c-numbers, anticipating 
a path integral treatment of the gauge fields. The time derivatives needed for Eq. (|2.10p can 
subsequently be taken according to the canonical equations of motion. 



(2.16) 



where the partial derivative operates on the background gauge fields. The commutator 
is readily evaluated with the help of equal-time anticommutators, and we also note that 
since Eq. (j2.10p includes a spatial integral over the currents, partial integrations are allowed. 
Adding together the two parts in Eq. ()2.16p then yields 



dT 
dt 



1 



M I 



(2.17) 



This can now be inserted into Eq. p.lOp . whereby the explicit factors of M duly cancel, since 
Mkin = M up to 0{T/M) thermal corrections which vanish in the heavy quark-mass limit: 



Elim 



M— >oo % 



1 

~00 



dt d'^x 



gE'4> - e^gE'e] (t, x) , \4>^gE'(i) - e^gE'o] (0, 0) 



(2.18) 

At this point the heavy quarks have become purely static; the ordering of the limits no longer 
matters, so we have set w — > inside the Fourier transform. 

Given that our derivation made no use of weak-coupling approximations, we believe that 
Eq. (j2.18p is free from (even finite) renormalization to all orders in perturbation theory, in 
the assumed regularization schemes with no additive mass renormalization and C2 equal to 
unity. This shows, in particular, that k is M- independent. 
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2.3. Euclidean correlator 



Eq. ()2.18p is a two-point function of gauge-invariant local operators; it therefore satisfies the 
standard KMS conditions which allow us to relate it to a Euclidean correlation function. In 
particular, let us define the Euclidean correlator 

Ge{t) ^ -f E Jim^^ U^^(^[(t>'^gEi(t>-e^gEie\{T,^) [cj^'^ gEict^-O^ gEie\{Q,Q)) . (2.19) 
i=i 

Hats have been left out because regular Euclidean path integral techniques apply for this 
object, and the minus sign accounts for the fact that a Euclidean electric field differs by a 
factor i from the Minkowskian one. The corresponding spectral function can be determined 
by inverting (for recent practical recipes see, e.g., refs. |14j^ the relation 



Ge{t) 



or analytically from 



GEi^n) 

p{uj) 



P(w) 

^ 







sinh^ 



dTe"^"^GE{T) , 

ImGEi^n — > —i[uJ + iO^]) 



(2.20) 



(2.21) 
(2.22) 



The momentum diffusion coefficient then follows from 

2T 

K = lim — piuj) . 

a;^0 id 

Note also that by making use of Eq. (|2.13|) , the susceptibility defined in Eq. 
the Euclidean theory be written as 



(0,0) 











= j\r j d3x( 




(r,x) 













(2.23) 



can m 



(2.24) 



In order to work out the contractions in Eqs. (j2.19p . (j2.24p . we need the heavy quark 
propagators within the Euclidean theory 



Ce = e\Dr + M)e + (l)\Dr - M)(j) + O 



1 

M 



(2.25) 



*The inversion leads to well-known systematic uncertainties, and we have nothing concrete to add on how 
to treat those. However, as will be demonstrated below, our spectral function is smoother at small uj than the 
ones in refs. QJ, which should somewhat ameliorate the problems in reconstructing the spectral function. 
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Making use of the equations of motion satisfied by the propagators, together with the proper 
boundary conditions, it can be shown that in the M — > oo hmit and for r > 0, 

= <5(3)(x-y)C/,^(r,0)e-^^^ 
= -5(3)(x-y)C/.^(/3,r)e(-/^)^, 



0,(r,x) 9*^(0, y) 
0„(O,x)^;(T,y) 
'Aq(t,x) 0^(0, y) 
'Aa(0,x) 4>1{t, y) 



-<5(3)(x-y)C/t^(r,0) 



(2.26) 
(2.27) 
(2.28) 
(2.29) 



where U is now a straight fundamental Wilson line in the Euclidean time direction. With 
these propagators, we obtain 



d^x 



^^gEi(t> - e^gEiO (r, x) gEicf^ - O^gEiO (0, 0) 



4(5(3) (0)e-^^^ (Re TV [[/(/3, r) gEi (r, 0) C/(r, 0) gE^ (0, 0)] 

,00 



(2.30) 



Similarly, the susceptibility x can be written as 



X 



00 



4(5(3) (0)e-'3^^ ^ dr (Re TV [;7(/3, r)C/(r, I 
4(5(3) (0)e-'^^^/3/ReTr [U{(3,0)]' 



In total, then. 



Ge{t) 



1 3 (Re TV 
3 ^ 



?7(Ar)ff^.(r,0) C/(t, 0)5^^(0,0) 



Re TV [[/(/?, 0)] 



(2.31) 



(2.32) 



and K can be obtained from the corresponding spectral function through Eq. (j2.23p . Note 
that the correlation function Ge{t) is positive (in a gauge with vanishing Aq, one can think 
of it as —drdaF{T — a)\a=o, where F is the correlation function of Ai). Eq. ()2.32p is our main 
result. A related formula in Minkowski signature was given in ref. [8j. 

It is appropriate to remark that the meaning of Eq. (j2.32p is unclear in the confinement 
phase of pure SU(A'c) gauge theory, where the expectation value of the Polyakov loop vanishes. 
In this situation, however, there would be a flux tube which drags the heavy quark in a way 
that is quite unlike diffusion, so it need not be surprising if the result for a diffusion coefficient 
were ill-defined. 



3. Perturbation theory 

The derivation of our main result, Eq. (j2.32p . made no use of the weak-coupling expansion, 
and is meant to be applicable everywhere in the deconfined phase, particularly at the phe- 
nomenologically interesting temperatures of a few hundred MeV. Nevertheless, to gain some 
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understanding on the general shape of the corresponding spectral function, we now go to very 
high temperatures, where the weak-coupling expansion is applicable. Our goal is to demon- 
strate explicitly that even in this regime, where spectral functions in general have more peaks 
and cusps than in a strongly-coupled regime, ours is relatively smooth. 

The leading-order (free theory) behaviors of the correlation function in Eq. (j2.32p and of 
the spectral function in Eq. ()2.22p are easily found: 



Ge{t) 
p{uj) 



+ 



1 



cos (vrrT) 

sin^(7rTr) ' 3sin'^(7rTr) 



Gvr 



(3.1) 
(3.2) 



where Cp = [N^ — l)/{2Nc)- This shows that at the free level the spectral function has no 
zero-frequency peak, in contrast to the spectral functions relevant for transport coefficients 
and vector current correlators (which have (5- function peaks at this order). Given that p{uj) 
in Eq. (13. 2p vanishes faster than oc uj, the diffusion constant k of Eq. (I2.23P is zero; we must 
work harder to find the leading non-trivial behavior at small frequency. 
At next-to-leading order, 0{g'^), the intercept k becomes non- vanishing [3]: 



g'CpT 2 A 2r 1 

— mjy In + - 

DTT V mD 2 



7e + 



C'(2) 



+ 



iVf In 2 



C(2) 2N, + N{ 



l + Oig) 



(3.3) 



where = g'^T'^{Nc/3 + Nf/6). As indicated, corrections to this expression start already 
at 0{g), and have in fact recently been determined [lOj . 

In order to learn how "easy" it is to extract the intercept k in practice, let us calculate 
more carefully the small-tij behavior of the spectral function in Eq. ()2.22p . We restrict, in 
the following, to frequencies at most of the order of the plasmon (or Debye) scale, uj < gT. 
Defining k{u)) to be the product on the right-hand side of Eq. (|2.23p . the difference [k{u)) — k] 
gets contributions only from soft momenta k ~ mj), and can be calculated at tree-level using 
Hard Thermal Loop propagators. Moreover, the Wilson lines in Eq. (|2.32p can be set to 
unity. Inserting the gluon propagator 



JK-{x-y) 



K 



(K) 



+ 



(3.4) 



where ^ is the gauge parameter, and carrying out the Fourier transform in Eq. (I2.2ip . we get 



g^Cp f d^k 



(2vr) 



2a;2 



+ 



ujl + k'^ + nT(u;„, k) u;2 + fe2 + Ilpiun, k) 



(3.5) 



where k = |k|. After the analytic continuation in Eq. p.22p . — > —i[uj + iO^], the self- 
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energies become (see, e.g.. ref. |12J ) 
UT{-i{uj + iO+),k) = 
nij(-i((j + iO+),k) = 



21 



CO 
l2 



In 



LU + iO^ + k 
uj + i0+ -k 



1-^ 

fc2 



cj , UJ + iO^ + k 

1 m ; 

2k LO + iO+ -k 



(3.6) 
(3.7) 



This leads to Landau cut contributions at k > u, and plasmon pole contributions at A; < w. 
Concretely, 

2g'^CFT Att 



k{u!) — K 



X vrmj-) X 



J U) 



dk r 



P - U)2 + i 



dk r 



— + — — lln ^"^'•^ 
k'^ 2k\ fc2 y fe-w 

1 



0(A; - (2)) X 4 X ^ 



2 / . \ 2 

if '1 



J_ 

2fc 



k2^l_u^l^k+uA ^(^>A (A:2 + l) 

2k k-Cj J \2k J 



|3(4 - (£;2)2 _ c2;2 



+ 



1 kU 



- kl 



UJ 



|3(4-(2;2) + l| 



fc|+l__^ln±tiB =0 



(3.8) 



2fe p LU — kr, 



where uj = w/mo and k = k/rriD. The four terms correspond to the transverse cut, electric 
cut, transverse pole, and electric pole, respectively. 

The outcome of a numerical evaluation of Eq. (j3.8|) is plotted in Fig. [U in units of the 
coefficient g^CpTm^/Q-K multiplying the logarithm in Eq. (|3.3|) . For uj < l/\/3, the result 
comes exclusively from the Landau cuts; for uj > plasmon poles contribute as well. 

For UJ ^ rriD, the result is dominated by the transverse pole, and extrapolates towards 

vr X 2a)2, the free theory result. 

The pattern in Fig. [T] illustrates an important point: even at weak coupling there is no trans- 
port peak around the origin; rather p{uj)/uj displays a relatively flat behavior at w < mD/\/3, 
with a significant rise only above the Debye scale. The only singularity is associated with 
the onset of the plasmon contributions; however this should be smoothed out in the full dy- 
namics. The amount of smoothening can be estimated from the (zero-momentum) plasmon 
damping rate calculated in ref. [15J, Fpi = Q.QAg'^NcT /2A'k. Already for = 0.05 this gives 
(for Nc = 3, A^f = 0...4) a width Fpi/mo > 0.2 comparable to that of the cusp; therefore we 
expect the true behavior to be completely regular. A more detailed study of the shape, in- 
cluding the effects of interactions in the small-u; regime and ultraviolet features in the large-w 
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Figure 1: A numerical evaluation of Eq. (|3.8p . in units of g^CFTm^j/GTr. The cusp is a feature of 
the weak-coupling expansion, as discussed in the text. 

regime, is deferred to a future publication ^16j . We also remark that the corresponding spec- 
tral functions computed for = 4 Super- Yang-Mills theory at infinite 't Hooft coupling show 
an analogous behavior, with the smooth infrared part ending in that case at w ~ T [8} [T7]. 

4. Correlator in lattice regularization 

Let us finally move to lattice regularization. In principle correlators of the type in Eq. (j2.32p 
can be measured with standard techniques on the lattice, in fact even at low temperatures 
where the signal is very small [18]. There is the problem, however, that the lattice electric 
fields require in general multiplicative renormalization factors (see, e.g., ref. [E]); these de- 
pend on the details of the discretization chosen, and it is also not clear how they could be 
determined on the non-perturbative leve|§. 

It appears, however, that the problem can at least be ameliorated if we choose a discretiza- 
tion of the electric fields inspired by lattice heavy quark effective theory (see, e.g., ref. |21j). 
The spatial components of the current (Eq. (j2.14p l could be thought of as 



J' 



I 



e\x + j)u]{x)e{x) - e\x)Uj{x)e{x + ]) -{§ — 



(4.1) 



2aM L 



^For recent progress with lattice magnetic fields, see ref. [20j . 
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where a is the lattice spacing; j is a unit vector in the j-direction; and Uj is a spatial link 
matrix. Discretizing also the time derivatives in Eq. (|2.10p and carrying out the contractions, 
we end up with a representation of Eq. (j2.32p which can best be represented graphically: 



(r-j) 



(L-n) + 



(L-n) 



(r-j)> 




Ge{t) 



6a^ ReTr( 



) 



(4.2) 



Here the direct lines within parentheses are link matrices; reading from the right, the long 
horizontal Wilson lines in the numerator have lengths t — a and (3 — t — a, if the sources are 
placed around xq = a/2 and xq = r + a/2, respectively; and the denominator stands for the 
trace of the Polyakov loop. It appears that Eq. ()4.2p should be less ultraviolet sensitive than 
the usual discretizations of the electric fields [181 US] • 

Continuing with the framework of the lattice heavy quark effective theory, the renormal- 
ization of Eq. (j4.2p can also be discussed in concrete terms, and be related to two separate 
issues. First of all, the linearly appearing mass parameter M in Eq. (|2.1ip is no longer the 
pole mass but requires additive renormalization; second, the coefficient C2 can differ from 
unity due to the absence of (Euclidean) Lorentz invariance. It appears that both of these 
issues could be addressed perturbatively and, in fact, even non-perturbatively [2T]. Since the 
explicit results depend on the particular lattice discretization chosen we do not, however, go 
into details here. 

5. Conclusions and outlook 

The main purpose of this paper has been to give a non-perturbative definition to the heavy 
quark momentum diffusion coefficient, k, allowing in principle for its lattice measurement. 
The basic definition is given in terms of a certain limit of the vector current correlation 
function, Eq. ()2.9|) . Making use of heavy quark effective theory, we have furthermore shown 
that the definition can be reduced to a much simpler purely gluonic correlator, given in 
Eq. K3% . with k then following from Eq. ({2:23]) . 

An important consequence of these relations is that they show that k, does not contain any 
logarithms of the heavy quark mass M. Our formulae could in principle also serve as the 
starting point for a first computation of a finite-temperature real-time quantity to relative 
accuracy Og, revealing in particular how the renormalization scale should be fixed. 

Moving to the non-perturbative level, we have also suggested a particular discretization 
of Eq. (j2.32p . given in Eq. (j4.2p . which could be free of significant renormalization factors. 
It remains to be tested in practice, however, how noisy the correlator is, and how fast the 
continuum limit can be approached. In addition, current practical recipes [14J related to the 



12 



inversion of Eq. (j2.20p suffer from uncontrolled systematic uncertainties which our method 
does not remove completely, although we hope that from the practical point of view they are 
less serious than in many other cases. 

Assuming that a non-perturbative value can be obtained for k, we can finally proceed to 
consider the thermalization rate of heavy quarks. A concrete and theoretically satisfactory 
meaning for a thermalization rate is provided by the heavy quark relaxation rate, or drag 
coefficient, denoted hy ijjj and defined around Eq. ()2.5p (the relation to thermalization follows 
from Eq. (jA.lh ). Employing the fluctuation-dissipation relation in Eq. (j2.8p . r]D can be 
estimated as r]£) ~ k/2MT, where M is the heavy quark pole mass and T is the temperature. 
Although this relation does have ambiguities related to the definition of the quark mass 
(a pole mass has inherent non-perturbative ambiguities at the level of several hundreds of 
MeV [22]; a treatment free of this problem can only be given in terms of non-perturbatively 
renormalized heavy quark effective theory [21]), such ambiguities should be subdominant 
compared with the large corrections related to infrared sensitive thermal physics, at least for 
the bottom quarks. These large thermal corrections are properly captured by our definition 
of K, so r]£) should lie in the right ballpark as well. We are therefore very much looking 
forward to the first numerical estimates of k. 
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Appendix A. Examples of dynamics leading to a transport peak 

In this appendix, we review briefly two arguments through which the Lorentzian form of the 
transport peak in Eq. (12. 4p can be established explicitly. 

Consider first non-relativistic quantum mechanics. Let us define Vi = Pi/M]^^^, where 
Pi is the momentum operator of the heavy quarks (i.e. the generator of translations in their 
Hilbert space). Suppose that we have, through some external source field, managed to prepare 
a non-equilibrium state where there is a heavy quark with a non-zero velocity. In thermal 
equilibrium, the average velocity must vanish, so we may expect the system to behave as 




non-eq — 



non-cq 





13 



Once t is so large that ('Vi(t)) non-eq ~ [{''^i)eq]^^'^ , Brownian motion sets in, and the system 
effectively equilibrates. In equilibrium we may define the correlator 

A,,it)^(^^{v,{t),vM}^ ■ (A.2) 

This is an even function of t and must vanish for t ^ oo; in fact, at least on certain time 
scales, it can be argued that it vanishes with the same exponent as the non-equilibrium 
correlator in Eq. (jA.ip (see, e.g., §118 of ref. [23]): 

A,i(t) ' ~'^A,,e-''^l*l , (A.3) 
where An is a constant. Taking a Fourier transform yields 

Aiiiu;)^ r dt e*-* Au (t) Au , (A.4) 

and making use of the general relation Aj.j(Lj) = [1 + 2?7,B(a;)] pii{<jj) (see, e.g., ref. [E]), where 
n-B{ijj) = l/[exp(/3cj) — 1] and pii{uj) is the spectral function, we arrive at 

to 2T uj'^ + rjfy 



This indeed agrees with the functional form of Eq. ()2.4p . 

Another example is given by classical Langevin dynamics (see also ref. [6]). Essentially, we 
replace {pi (t)) non-eq Pi{i)i and assume the dynamics to be contained in 

Pi{t) = -VDMt)+Ut), (A.6) 
HmC.it'))) = K,,5,,5{t-t'), {{m))=0, (A.7) 

with ^ a Gaussian stochastic noise field, and ((...)) denoting an average over the noise. For a 
heavy particle the Gaussian nature follows from the central limit theorem and the slow time 
scale of its dynamics, while the auto-correlator Kd = /^dt ((?j(i)'^j(0))) can be chosen such 
as to match that of the underlying theory, Eq. (j2.5p . It is easy to verify that within this 
dynamics, for a distribution with density Tx^^ of heavy quarks, the equilibrium correlator is 
exactly Eq. p7 
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